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Abstract. I review the present status of two related models addressing scenarios in which formation of 
heavy quarkonium states in high energy heavy ion collisions proceed via "off- diagonal" combinations of a 
quark and an antiquark. The physical process involved belongs to a general class of quark "recombination" , 
although technically the recombining quarks here were never previously bound in a quarkonium state. 
Features of these processes relevant as a signature of color deconfinement are discussed. 
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1 Introduction 

The original idea of Matsui and Satz [J predicted a sup- 
pression of J/ip produced in heavy ion collisions as a re- 
sult of the expected screening of the color force above the 
deconfinement phase transition. The prediction of sup- 
pression follows from the expectation that the eventual 
hadronization of the dcconfined charm quarks is preferen- 
tially with light up and down quarks, since generally only 
one cc pair is produced in a given collision. 

Several years ago, it was pointed out that the sup- 
pression scenario could be altered in nuclear collisions at 
collider energies PlEI- At sufficiently high energy, mul- 
tiple pairs of heavy quarks will be produced in a single 
nucleus-nucleus collision. Then it may be possible for a 
given heavy quark to form a heavy quarkonium hadron 
by combining with a heavy antiquark which originated 
from a different initial production process. I will refer to 
such combinations as "off-diagonal" pairs. The probabil- 
ity to form heavy quarkonium will of course depend on 
the physics of the interaction and also the nature of the 
medium in which it occurs. However, one can predict a 
few simple properties of this formation process based on 
general considerations. 

We consider scenarios in which the formation of J/ip 
is allowed to proceed through any combination consist- 
ing of one of the Nc charm quarks with one of the iVg 
anticharm quarks which result from the initial produc- 
tion of Ncc pairs in a central heavy ion collision. For a 
given charm quark, one expects then that the probability 
V to form a J/^p is just proportional to the number of 
available anticharm quarks relative to the number of light 
antiquarks. 



number of available charm quarks Nc to obtain the total 
number of J ftp expected in a given event. 



(1) 



where we normalize the number of light antiquarks by the 
number of produced charged hadrons. Since this probabil- 
ity is generally very small, one can simply multiply by the 



Nj/^ oc Nc, 



'/Nch, 



(2) 



where the use of the initial values Ncc = Nc = Nc is again 
justified by the relatively small number of J/ip formed. For 
an ensemble of events, we calculate the average number of 
J/V' per event from the average value of Ncs^ , and neglect 
fluctuations in Nch- 

< Nj/^ A(< Ncc > +1) < Ncc > /Nch, (3) 

where we place all dynamical dependence in the parameter 
A. The resulting quadratic dependence on < Ncc > Pro- 
vides a unique signature which must at some high energy 
become dominant over production via a superposition of 
independent diagonal cc pairs. 

Initial estimates of Ncc ~ 10 for central Au-Au colli- 
sions at RHIC used extrapolations of cross section mea- 
surements at lower energy More recently there are 
measurements at RHIC based on high-transverse momen- 
tum electrons by PHENIX and also reconstructed D-mesons 
by STAR which imply larger numbers. Central values of 
these measurements lead to Ncc « 20 (PHENIX) 5 or 40 
(STAR)[B), with relatively large experimental uncertain- 
ties which leave the two measurements consistent. In the 
following estimates for J/ ip we explore this entire range of 
initial Nrr- 



2 Statistical Hadronization 

This model was motivated by the success of predictions for 
the relative abundances of light hadrons produced in high 
energy heavy ion interactions in terms of a hadron gas in 
chemical and thermal equilibrium. Such fits, however, un- 
derpredict the abundances of hadrons containing charm 
quarks. This can be understood in terms of the long time 
scales required to approach chemical equilibrium for heavy 
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quarks, starting from the large number of charm quarks 
produced via hard processes during the initial stages of 
the collision. The original formulation 2 of the statis- 
tical hadronization model for hadrons containing charm 
quarks assumes that at hadronization the charm quarks 
are distributed into hadrons according to chemical equi- 
librium, but adjusted by a factor 7^ which accounts for 
oversaturation of charm. One power of this factor multi- 
plies a given thermal hadron population for each charm or 
anticharm quark contained in the hadron. Thus the rela- 
tive abundance of J /tp to that of D mesons, for example, 
will be enhanced in this model. The enhancement factor 
is determined by conservation of charm, again using the 
time scale argument to justify neglecting pair production 
or annihilation before hadronization. 



thermal model only in this region.) However, the charm 
production cross sections used in these calculations were 
different (390 /ib vs. 650 /ib), which would imply a differ- 
ence in predictions of almost a factor of 3, all other effects 
being equal. Although the thermal parameters appear to 
be compatible, the extraction of a rapidity-density volume 
parameter evidently is different in these two approaches. 



2.2 Centrality dependence 

It is now conventional in heavy ion collisions to parameter- 
ize the centrality of the collision in terms of the number 
of nucleon participants, Npart- The point-like process in 
which cc pairs are produced then leads to a 4/3 power law 
behavior of Nrr. 



^cc — (^'^ c-^ opencharm ~t~ Tc ^hiddencharnij (4) 

where Nopencharm and Nhiddencharm are calculated in the 
thermal equilibrium model grand canonical ensemble. (For 
peripheral collisions the total particle numbers are not 
sufficiently large and one must calculate in the canoni- 
cal ensemble.) It was first shown numerically in Ref. 0, 
and later motivated in Ref. |S] that the canonical correc- 
tion effect is equivalent to directly using the grand canon- 
ical value for 7c in the ensemble average according to Eq. 
13 The hidden charm term is negligible for all cases of 
present interest, and one finds that 7c is directly propor- 
tional to Ncc- Then one can express the number of J/^p at 
hadronization as 



(5) 



where n and V are the number density and volume ap- 
propriate to the relevant hadronization region. Insertion 
of the expression for 7c from Eq. ^ leads to an expression 
which has the form expected in Eq. |31 with 



in, 



ij^open)'^ 



(6) 



(Note that the factor of Nch appears due to replacing the 
one remaining power of V by the ratio of total number to 
density for charged particles.) 



2.1 Comparison with RHIC measurements 

There has been one measurement at RHIC by PHENIX 
for J/ -0 production in Au- Au interactions at 200 GeV [H] . 
The data was analyzed in three centrality regions, but 
due to limited statistics the uncertainties were quite large. 
Also, the most central data leads to only an upper limit. 
Two separate groups |10llH] have applied the statistical 
hadronization model in this case. Both have found general 
agreement with the data, which involves rapidity densi- 
ties at y = rather than total yields. (Only the two more 
central data points can be used, since measurements of 
the relative yields of J /ip and ip' are consistent with the 
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^part 

2A 



4/3 



(7) 



which is normalized by the maximum number A^cc(O) at 
impact parameter b=0 where N, 



part 



2k. 



One also requires the centrality dependence of Nch: for 
which we parameterize Nch = o,{NpartY~^'^ , where a is a 
normalization factor and Z\, which depends on the produc- 
tion process for charged particles, will be varied. Then the 
centrality dependence of Nj/^, is determined by substitu- 
tion into Eq. One generally normalizes the experimen- 
tal yield of J/ "tp by either the number of binary collisions 
(equivalently TVce), or the number of nucleon participants 
Npart- For the first normalization choice, one obtains 



^ (iNpart)^ 



Nr.. 



-A 



(N, 



part ) 



(8) 



This combination of two power-law terms in Npart 
which differ by | is obviously due to the combination of 

quadratic and linear terms in Ncc for Nj f^, . ^j^^ will have 

a minimum (for A < i) at Npart = N„iin = 2 A ^ i 
The sharpness of the minimum can be characterized by the 
ratio 



-L^ AT V ^ 



R = 



)/^^{2A), with the result 



4(1 -Zi) 



(zi-l/4) 



3(iVcc(0) + 1) 



NcciO) 



l + A 



(9) 



These features are shown in Fig. ^ for a range of A 
and iVcc(O). Aside from the curves for Z\ = 1 which are 
constant for large Npart, all of the minimum points are 
at relatively low values of Npart, in the region where the 
statistical hadronization cannot be applied. 

Figure 121 shows the corresponding behavior for the ra- 



tio 



N 



-, which can be obtained by making the substitu- 
tion A ^ A — ^. The same general behavior is seen, but 
the sharpness of the approach to the minimum value is 
enhanced, especially for the largest values of 7Vcc(0). The 
real test of the predicted centrality behavior requires data 
on both the magnitude of the initial charm production 
and the centrality dependence of Nch- 
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will exist in an environment at temperatures well above 
the deconfinement transition jlllll2| . The dominant for- 
mation process in this scenario involves the capture of a 
quark and antiquark in a relative color octet state into the 
color singlet J/ tp with the emission of a color octet gluon. 
This reaction is just the inverse of the primary J/tp disso- 
ciation process via collisions with deconfined gluons 'IHj. 
One can then follow the time evolution of charm quark and 
charmonium numbers in a region of color deconfinement 
according to a Boltzmann equation in which the formation 
and dissociation reactions compete. 
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Fig. 1. Centrality dependence of binary-scaled J/^ formation 
via statistical hadronization. 



dN 



dr 



XfNcNc[V{t)]-^ ^XuNjf^ Pg. 



(10) 



with Pg the number density of gluons. The reactivity A is 
the reaction rate (cruroi} averaged over the momentum dis- 
tribution of the initial participants, i.e. c and c for Xp and 
J and g for Ad. The gluon density is determined by the 
equilibrium value in the QGP at each temperature, and 
the volume must be modeled according to the expansion 
and cooling profiles of the interaction region. 

This equation has an analytic solution in the case where 
the total number of formed J /ip is much smaller than the 
initial number of N^c- 
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Fig. 2. Centrality dependence of participant-scaled J/^ji for- 
mation via statistical hadronization. 



3 Kinetic formation 

The kinetic model|3] describes a scenario in which the mo- 
bility of initially-produced charm quarks in a space-time 
region of color deconfinement allow formation of quarko- 
nium via "off-diagonal" combinations of quark and an- 
tiquark. The motivation for such a scenario in the case 
of J /ij} formation has received support from recent lattice 
calculations of spectral functions. These indicate that J /ip 



Nj/^Tf) = eiTf)[Nj/^iTo)+N!, r Af [V{t) e(r)]-i dr], 

J To 

(11) 

where Tf and tq are the final and initial times. The func- 
tion e(T/) — e would be the suppression fac- 
tor in this scenario if the formation mechanism were ne- 
glected. 

One can readily see that Nj/^ obeys the general prop- 
erties present in Eq. |31 However, the factor equivalent to 
system volume V (in the statistical hadronization model) 
is time dependent and modified by a combination of fac- 
tors involving the interaction rates. Thus the centrality 
behavior will depend on additional parameters. The ini- 
tial calculations [7] used the ratio of nucleon participants 
to participant density to define a transverse area which de- 
fines the boundary of the region of color deconfinement. 
This is supplemented by longitudinal expansion starting 
at an initial time tq — 0.5 fm (Transverse expansion was 
initially neglected, but has been included in subsequent 
calculations JH.) The expansion was taken to be isen- 
tropic, which determines the time evolution behavior of 
the temperature. The initial value Tq is taken as a pa- 
rameter, and the final Tf is fixed at the hadronization 
point. The reactivities Xp and A^i require specification of 
cross sections. For ao we use the "OPE-inspired" model 
of gluon dissociation of deeply-bound heavy quarkonium 
[T^ Pj3|, which is related via detailed balance to the cor- 
responding ap- These cross sections are shown in Fig. |31 
One sees that they are peaked at low energy, and that 
Up > gd due to the large binding energy. 

Fig. 21 shows the generic behavior of the J/ ip popula- 
tion resulting from a numerical solution of Equation 1101 
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t (arbitrary units) 

Fig. 4. Time development of J/V' with formation and dissoci- 
ation rates. 



systematic uncertainties. The lines shown are calculations 
in the kinetic model with a range of principal parame- 
ters, including Ncs (b=0), Tq, transverse expansion veloc- 
ity (vtr), and in one case an initial population fraction (x) 
of J /tp. The range of these parameters was chosen to ex- 
hibit what constraints are placed by this initial data. All 
of the calculations used the same charm-quark distribu- 
tion, which was taken from a LO pQCD calculation 
One sees that there is a substantial range of parameters 
allowed by this data, but that the increase in statistics an- 
ticipated in run 3 will allow a much more stringent limit 
for the acceptable (if any) region of parameter space in 
the kinetic model. 

The kinetic model also makes predictions for the mo- 
mentum space distribution of formed J/^. For this pur- 
pose we require the differential cross sections related to 
the (Tpjj. These are obtained via an adaptation of the 
corresponding expressions for photodissociation of atomic 
bound states. One can then express the time-integrated 
formation rate in terms of a sum over all cc pairs, each 
weighted by differential formation probabilities. 



One can see that the final population is in fact determined 
by the time integral of the difference between formation 
and dissociation rates, shown as dashed lines. The magni- 
tudes are determined by the parameter Tq, which controls 
the magnitude and time dependence of the gluon density, 
and also the total lifetime. It is important to note that for 
parameter values in a range consistent with expectations, 
the expansion rate of the color deconfinement volume and 
the decrease of gluon density with time prevent the sys- 
tem from reaching an equilibrium population within this 
lifetime. 

Fig. |S1 shows the PHENIX data for J/ip production in 
Au- Au collisions at 200 Gev 'S" . The most central bin yield 
only allowed an upperlimit (hatched horizontal lines), while 
two less central bins yielded absolute values plus addi- 
tional one-sigma upper limits both from statistical and 



(12) 

Note that the formation magnitude exhibits the explicit 
quadratic dependence on total charm, normalized by a 
prefactor which is proportional to the inverse of the system 
volume. 

We first look at the rapidity spectra of cc pairs, shown 
in Fig. |B1 and compare with the measured J/^p distri- 
bution in pp interactions |17j . Normalized spectra are 
used throughout, so that the results are independent of 
the prefactors. One sees that the data are consistent with 
the distribution of unbiased diagonal pairs only, which 
is what one would expect for pp interactions. The dis- 
tribution of all pairs (also unbiased) is somewhat nar- 



R. L. Thews: Quarkonium formation in statistical and kinetic models 



5 




rower than the data. Also shown are all pairs for which 
each is biased by the total formation probability appro- 
priate for the given pair energy. It is very close to the 
next curve, which takes into account formation from all 
pairs, using exact J /ip kinematics and the full differen- 
tial dependence. One sees that both of these curves are 
substantially narrower than the pp data. Thus the ki- 
netic model predicts that the rapidity distribution of J /ip 
formed by off-diagonal pairs (only possible in nucleus- 
nucleus collisions which lead to color dcconfinement) will 
be substantially narrower than J /tp produced in pp in- 
teractions at the same energy. Fig. [3 shows the trans- 
verse momentum spectra of unbiased diagonal cc pairs, 
along with the PHENIX data ^3 for J/i/; production 
in pp interactions at 200 GeV. The set of curves result 
from augmenting the quark initial momenta with a trans- 
verse momentum "kick" to simulate confinement and ini- 
tial state effects. The pp data restricts the magnitude of 
this kick, parameterized by a Gaussian distribution, to 
lie within the range < fcf >pp= 0.5 ± Q.lGeV"^. To ex- 
tend this to formation in Au-Au collisions, we must ex- 
tract the appropriate kt for initial state effects in the nu- 
cleus. We use PHENIX data for J /ip in d-Au collisions 
|18j , which shows that the pt spectra are broadened rela- 
tive to that in pp interactions. This results in an estimate 
for < k'l >Au-Au= 1-3 ± O.SGeF^, where the uncertainty 
is set by the rapidity variation of the J/ip pt broadening. 
This range of values was utilized in the formation calcula- 
tions in Au-Au interactions. The predicted rapidity spec- 
tra are found to be essentially independent of the magni- 
tude of the initial charm quark kt kick, so that the narrow- 
est of the curves in Fig. will serve as the kinetic model 
prediction for J/tp formed in an Au-Au collision. Fig. |S1 
shows the predicted transverse momentum spectra of J/ip 
at RHIC which would result from the formation mecha- 
nism, using the entire allowed range of < k^ >au-Au- For 
comparison we show the distribution of diagonal unbiased 
cc pairs with the central value in the allowed range of ini- 
tial < fcj >Au~Au, which should be relevant if all of the 




Pt (GeV) 



Fig. 8. Prediction for J/ip pt distribution from formation pro- 
cess in Au-Au collisions at RHIC. 

J/tp were produced directly from the initial cc pairs. Of 
course, both of these distributions would be modified by 
the competing dissociation process during the expansion 
phase, but one would anticipate a similar effect on each 
which would preserve the relative comparison. (A sample 
suppression factor applied to these curves actually shows 
very little change in the shape of the normalized spectra.) 



4 Summary 

One expects on general grounds that heavy quarkonium 
production in high energy heavy ion collisions must con- 
tain a component which is formed either during a period 
of color dcconfinement or at the hadronization point. The 
magnitude of this formation will increase quadratically 
with the total amount of charm initially produced via 
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nucleon-nucleon interactions. Both the statistical hadroniza- 
tion model and the kinetic formation model exhibit this 
property. The absolute magnitude is somewhat model- 
dependent, and initial RHIC data for J /ip can be accom- 
modated. The rapidity and transverse momentum spec- 
tra may be decisive in determining whether or not this 
formation makes a significant contribution. The kinetic 
formation model predicts a narrowing of the J /ip rapid- 
ity distribution (compared with that in pp collisions) , and 
also a narrowing of the transverse momentum distribution 
(compared with an extrapolation of behavior measured in 
pp and d-Au collisions). In principle, both of these forma- 
tion mechanisms can coexist, so that the upcoming Au-Au 
data may reveal a two-component structure. 

This work was supported by U. S. Department of Energy Grant 
DE-FG02-04ER41318. 
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